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Abstract 

This paper investigates the performance of optimal single stream beamforming schemes in multiple-input 
multiple-output (MIMO) dual-hop amplify-and-forward (AF) systems. Assuming channel state information is not 
available at the source and relay, the optimal transmit and receive beamforming vectors are computed at the 
destination, and the transmit beamforming vector is sent to the transmitter via a dedicated feedback link. Then, a 
set of new closed-form expressions for the statistical properties of the maximum eigenvalue of the resultant channel 
is derived, i.e., the cumulative density function (cdf), probability density function (pdf) and general moments, as well 
as the first order asymptotic expansion and asymptotic large dimension approximations. These analytical expressions 
are then applied to study three important performance metrics of the system, i.e., outage probability, average symbol 
error rate and ergodic capacity. In addition, more detailed treatments are provided for some important special cases, 
e.g., when the number of antennas at one of the nodes is one or large, simple and insightful expressions for the 
key parameters such as diversity order and array gain of the system are derived. With the analytical results, the 
joint impact of source, relay and destination antenna numbers on the system performance is addressed, and the 
performance of optimal beamforming schemes and orthogonal space-time block-coding (OSTBC) schemes are 
compared. Results reveal that the number of antennas at the relay has a great impact on how the numbers of 
antennas at the source and destination contribute to the system performance, and optimal beamforming not only 
achieves the same maximum diversity order as OSTBC, but also provides significant power gains over OSTBC. 
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I. Introduction 

Multiple-input multiple-output (MIMO) antenna technology can be used to substantially increase the 
spectral efficiency or significantly improve the reliability of the communication systems [1,2], and has 
become one of the key enabling technologies for future wireless communication systems. Parallel to the 
development of MIMO technology, dual-hop relaying technology, where an intermediate relay node helps 
forward the source signal to the intended destination node, is another efficient approach to considerably 
improve the performance and extend the coverage of the communication systems [3, 4]. The huge potential 
of achieving both the MIMO and relaying benefits has resulted in enormous interest in the MIMO relaying 
systems [5,6]. 

Due to the low complexity of amplify-and-forward (AF) relaying protocol, it is of great importance to 
understand the performance of MIMO dual-hop AF systems. Prior works on MIMO dual-hop AF systems 
have largely focused on the asymptotic characterization of the information theoretical measures such as 
capacity scaling behavior and diversity-multiplexing tradeoff (DMT). For instance, in [7], the impact of 
the number of antennas at the relay on the asymptotic capacity was studied, while [8, 9] investigated the 
asymptotic capacity of multi-hop AF relay systems using tools from large-dimensional random matrix 
theory. The DMT of MIMO dual-hop AF systems was established in [10], and was later extended to the 
multi-level AF systems [11]. 

A few results on MIMO dual-hop AF systems in the finite regime have also emerged recently. In 
particular, an exact ergodic capacity analysis was presented in [12], while in [13, 14], the symbol error 
rate (SER) performance of orthogonal space-time block-coding (OSTBC) schemes in such system was 
addressed using a moment generating function based approach. The performance of optimal beamforming 
schemes has also been investigated in various multiple antenna dual-hop AF systems [15-18]. However, 
the major limitation of all these prior works [15-18] is that one or more nodes is restricted to have only 
one antenna. In addition, all these works consider the channel state information (CSI) assisted relaying 
protocol, which however requires the relay having access to the CSI. With the implementation of multiple 
antennas, acquiring the CSI at the relay entails higher implementation complexity and cost, hence, it may 
not be suitable for certain applications with low cost budgets. On the other hand, the fixed-gain relaying 
protocol, where the relay simply forwards a scaled version of the received signal, becomes more appealing 
for such applications. Despite its practical importance, the performance of optimal beamforming schemes 
in general MIMO dual-hop AF systems with fixed-gain relaying protocol remains unknown. 

Motivated by this, we adopt the same system model as in [8, 12, 14] and study the analytical performance 
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of optimal beamforming schemes in the general MIMO dual-hop AF systems. It is worth pointing out 
that, with multiple antennas at all the three nodes, it is possible to send multiple independent date 
streams at the same time, namely, multi-stream beamforming systems (See, i.e., [20,21] and references 
therein). In general, the multi-stream beamforming system provides higher data rate, but experiences 
worse error performance compared with the single stream beamforming system. In the current paper, we 
focus exclusively on the single stream beamforming system, and leave the analysis of the multi-stream 
beamforming system for future works. In addition, we assume that the CSI is not available at the source and 
relay. Instead, it is estimated at the destination via the methods proposed in [22-24]. Hence, the optimal 
receive and transmit beamforming vectors are computed at the destination, and the transmit beamforming 
vector is then forwarded to the transmitter through a dedicated feedback link. 

The main contribution of the paper is the derivation of a set of new statistical results pertaining to 
the maximum eigenvalue of the underlying MIMO dual-hop AF channel matrix. In particular, we present 
exact closed-form expressions for the cumulative distribution function (cdf), probability density function 
(pdf), and general moments, as well as asymptotic first order expansion of the cdf and pdf of the maximum 
eigenvalue. In addition, the statistical behavior of the maximum eigenvalue is also characterized in the 
asymptotic large antenna regime. We then apply these statistical results to study the performance of optimal 
beamfomring systems in terms of three key performance metrics, i.e., outage probability, average SER 
and ergodic capacity. Furthermore, we provide detailed treatments for certain important special cases, 
and obtain key performance indicators such as diversity order, array gain and high signal-to-noise ratio 
(SNR) slope and power offset of the system. Finally, we compare the results of optimal beamforming 
with OSTBC schemes, and show that optimal beamforming schemes achieve the same maximum diversity 
order as OSTBC schemes, but have significant power gains over OSTBC schemes. 

The remainder of the paper is organized as follows. Section II introduces the system model. Section III 
presents a host of new statistical results of the maximum eigenvalue, while Section IV gives a detailed 
investigation on the performance of optimal beamforming schemes in terms of outage probability, average 
SER and ergodic capacity. Finally, Section V concludes the paper. 

Throughout the paper, the following notation is adopted. Vectors are represented as columns and are 
denoted in lower-case bold-face, while matrices are represented in upper-case bold-face. The superscript 
indicates the matrix conjugate-transpose operations, respectively. We use T(x) to denote the gamma 
function, det(-) to denote the matrix determinant, || • || to denote the Frobenius norm, E{-} to denote 
the expectation operation, I r to denote an r x r identity matrix, and CAf axb to denote an a x b matrix 
with entries being identically and independently distributed (i.i.d.) zero mean circular symmetric complex 
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Gaussian (ZMCSCG) random variables with unit variance. 

II. System Model 

We employ the same AF MIMO dual-hop system model as in [8, 12, 14]. Suppose that there are n s 
antennas at the source terminal, n r antennas at the relay and antennas at the destination terminal, 
which we represent this by the 3-tuple (n s ,n r ,n d ). All terminals operate in half-duplex mode, and as 
such communication occurs from source to relay and from relay to destination in two separate time slots. 
The end-to-end input-output relation of this channel is then given by 

y = H 2 FHiS + H 2 Fn nr + n nd , (1) 

where s is the transmit symbol vector, n, tr and n nd are the noise vectors at the relay and destination 



terminals, respectively, F = J - ^ +p ^ - nr ( a corresponds to the overall power gain of the relay terminal) 
is the forwarding matrix at the relay terminal which simply forwards scaled versions of its received 
signals!, and Hi E CJ\f nrXns and H 2 E CJ\f ndXnr are the MIMO channel matrices from the source to 
relay and from the relay to destination, respectively. The additive noise at the relay and destination are 
assumed to be white in both space and time and are modeled as ZMCSCG with unit variance, i.e., 
E { n n r i4 r } = l rir and E{n„ d nt d } = I„ d . 

For the transmit beamforming scheme, we have s = w t x where w t is the beamforming vector with unit 
norm, i.e., ||w t || 2 = 1, and x is the transmit symbol satisfying E{|x| 2 } = p. At the destination node, the 
receiver combines the signal by multiplying the received signal vector with w r . Hence, the post-processed 
signal at the receiver node can be expressed as 

wjy = y/awlH 2 Hiw t x + A/aw r t H 2 n nr + wjn nd , (2) 

where a = — -r—^- As a result, the end-to-end SNR is given by 

n r (l+p) ' b J 

wjH 2 HiWtw t t H{H 2 w r ap 

7 = — 77 1 \ • (3) 

w| ^aH 2 H 2 + l nd j w r 
It can be easily observed that the optimal combining vector w r in maximizing the SNR is 

w r = (aH 2 H^ + I nd ) _1 H 2 H lW4 . (4) 

'it is worth pointing out that, with the implementation of multiple antennas, the variance of the received power per antenna at the relay 
is reduced. 



Substituting © into ©, the corresponding maximum SNR is given by 

7 = apw\H{Hl (aH 2 H 2 + I„ d ) H 2 H lW< . (5) 

Therefore, the optimal transmit beamforming vector corresponds to the eigenvector associated with the 
maximum eigenvalue of H{H 2 ^ctH 2 H 2 + I„ d j H 2 H l5 giving the resultant SNR as 

7 = apA max , (6) 

where A max is the maximum eigenvalue of the random matrix H|H 2 ^aH 2 H 2 + I„ d j H 2 Hi. 

III. Statistical Properties of A max 

This section presents a set of new statistical results of A max , including exact expressions for the cdf, pdf, 
and general moments of the maximum eigenvalue A max , as well as the asymptotic first order expansions. 
In addition, the statistical behaviors of A max are also investigated in the asymptotic large antenna regimes. 
For notation convenience, we find it useful to define q = min(nd, n r ), p = max(nd,n r ), s = min(n s ,g), 
t = max(n s , q), m = min(n s ,p), n = max(n s ,p), 6(i,j) = 2q+p — i—j — s, and r(i,j) = s + i+j — q — 2. 

A. Exact Expressions for the cdf 

The following theorem gives the cdf of A max and will be useful in the analysis of the outage and SER 
performance of optimal beamforming scheme in MIMO dual-hop AF systems. 

Theorem 1: Let Hx £ Qn r y.n s an( j jj 2 ^ Qn d xn r ^ an( j a b em g a non-negative constant, the cdf of the 

maximum eigenvalue A max of HjH 2 ^aH 2 H 2 + l nd \ H 2 Hi is given by 

(-l)^-^det(<fr(x)) 

Xm ° x{) n?=i r (?-^ + i)r(p-z + i)' {) 

where is a q x q matrix with entries 

Effo { T{ ] j) )a l T(8(i,j) + 1 + 1)- e- ax B(x) for % > q - s. 



where 

t-i k r(i,j)+k . . . . N 
fc=0 1=0 

where K v (x) is the modified Bessel function of the second kind [25]. 

Proof: See Appendix HI □ 



t-i k r(i,j)+k , . . .. . 

B{x)=2\— }^ I V ' J; ja'x a lf flW)+ ,_ fe+1 (2 > /i), (9) 

u—n ' l—n \ / 



Fig.fflplots the cdf curves of A max with different system configurations. Results in this figure demonstrate 
that the analytical results are in perfect agreement with the Monte Carlo simulation results, which confirms 
the correctness of the analytical expressions. 

The following corollaries present two special cases for which the cdf expression in © can be reduced 
to a simple form. The results will be useful to study the outage probability of the systems where one of 
the nodes is equipped with only a single antenna. 

Corollary 1: When n s = 1, © reduces to 

n?=ir(<z-« + i)r(p-i + i 



^A max (x) = } - ; - ; 1N (10) 



where r] = p + q + l — j, and Dj is the (q, j)-th cofactor of a q x q matrix whose (i,j)-th entry is given 

by (_l)*-l-< (^"^oTCp + 1 + 

Proof: When n s = 1, we have s = 1 and i = g. As a consequence, the desired result follows by 
substituting s = 1 into © and then applying Laplace's expansion along the last row of the determinants, 
which completes the proof. □ 
Corollary 2: When q = 1, © reduces to 

= i - w E F E (*V^;W(2V5). (i i) 

fc=o i=0 ^ ' 

Proof: The proof is straightforward and thus omitted. □ 

5. Exact Expressions for the pdf 

The following theorem presents the pdf of the maximum eigenvalue A max . This will be used to compute 
the ergodic capacity of optimal beamforming scheme in MIMO dual-hop AF systems. 

Theorem 2: Let Hi e C' vx " s and H 2 G C ndXnr , and a being a non-negative constant, the pdf of the 
maximum eigenvalue A max of H^H^ ^aH 2 Hj> + I nd j H 2 Hi is given by 



Yl 9 i=l V(q-i + l)V(p-i + l) > ^ 



where is a g x g matrix with entries 



[<J>(x)L i for % ^ /, 

*i(»)]ii= { (13) 
2e— EU, A Effo J)+fc ( T(< f + *)o'y(x) for z = I. 



where 

6(i,j) + l + k-l 6(i,j) + l + k 

g(x) = (ax-k)x 2 K 6{i j )+ i_ k+1 (2y/x) + x 2 K e{i j )+l _ k (2y/x). (14) 
Proof: See Appendix flU □ 
Similar to our cdf analysis before, the following corollaries present two special cases for which the 
pdf expression in ([T2T ) is simplified. The results will be used to study the ergodic capacity of the systems 
where one of the nodes is equipped with only a single antenna. 
Corollary 3: When n s — 1, (PT2l) reduces to 

/AmJx) ~ IlLir(g-* + i)r(p-i + i) ' (15) 

where 77 and Dj are defined in Corollary [T] 

Proof: When n s = 1, we have s — 1 and t = g. Hence, the desired result is obtained by substituting 
s = 1 into (fT2l) and using Laplace's expansion on the last row of the determinants. □ 
Corollary 4: When q — 1, (fT2l) reduces to 



/A max (s) 



r(p) "f||(|)« H (H) , -V<« + ' 4! V,p4 (16) 

Proof: The proof is straightforward and thus omitted. □ 



C. Asymptotic Expansions for the cdf and pdf 

Here, we present some asymptotic expansions for the cdf and pdf of A max . Although obtaining such 
expansions for general systems with arbitrary n s , n r and rid is a challenging task, simple expansions are 
possible for the important special cases where one of the nodes has only one antenna, as we show in the 
following. These results will be used for investigating the asymptotic outage probability and for deriving 
the diversity order and array gain of the systems. 

Theorem 3: When n s = 1, the cdf and pdf of A max have the asymptotic expansions: 

FxUx) = -x q + o(x q ), (17) 
Q 

/A max (x) = ^- 1 + o(^- 1 ), (18) 

where 

detQg) 

Vl ~m m=ir(p-* + i)r(g-» + i)' ( } 
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where ^ is a q x g matrix with entries 

[*L'= 9 E ( 2? ~*~ J + 1 )a*r(p-g + * + i + j-2). (20) 

fc=0 

Proo/: See Appendix [Till □ 
Theorem 4: When q = 1, the cdf and pdf of A max have the asymptotic expansions: 

^A max (*) = -* m + o(* m ), (21) 
m 

fx^(x)=V2X m ~ 1 + o(x m - 1 ), (22) 

where 



( Ego (?)aT( P -» 3 +0 

r(p)r(n.) 

-c-lnx+ggi (° s )»T(i) 

r(p)r(n.) 

T(n a -p) 



if p > n s , 

if p = n s , (23) 
if p < n s . 



r(p)r(n.) 

Proo/: See Appendix |TVl □ 



D. Exact Expressions for the Moments 

In this section, we present exact expressions for the moments of A max for the important special cases 
corresponding to the systems where one of the nodes is equipped with a single antenna. These results 
will be used for investigating the ergodic capacity of the system. 

Theorem 5: When n s = 1, the m-th moment of A max is given by 



vU\ W (-l^lXm + l) « 

U maxj 1 ~ UU r(g - i + i)r(p - < + 1) ^ j 



E C Z 1 ) a ^ mr ^ + m ) (^ + ™)^ ( m+1 ' 1_ ^'~) + a_lf/ ( m+1 ' 2_7 7'~)) J • (24) 
Proof: Starting from the definition of the moment, the desired result can be obtained with the help 

of the following integration relationship 

x»e- mx K v (2^)dx = ^ I 2 >-V U - V - + 1, 1 - v, E- . (25) 

2m M 2+ 1 V 2 ml 



□ 



Theorem 6: When q = 1, the m-th moment of A max is given by 



E{(A max ) m } = E M F(P+ : + m) ^ (n s + m, 1 - P + n s - i, ±V (26) 

u max; / r(p)r(n s ) ^ V ^ / a na+m ~ l \ a J 

Proof: The proof is similar to that of Theorem [5] □ 



E. Asymptotic Approximations for the cdf 

In this section, we present the asymptotic approximation for the cdf of A max in the case where the number 
of antennas at one of the nodes becomes large. These results are useful when studying the average SER 
and ergodic capacity of the system. 

Theorem 7: When the number of antennas at one of the nodes becomes large, A max is statistically 
equivalent to 

,KAm - 1 ' 1 if n. -> oc. 



Amax 



aA max i+l s ' 



A m ax2 if n d 00, (27) 



-A maX 3 if n r -)> oo 



an r +l ' 

where A max i, A maX 2 and A maX 3 are the maximum eigenvalues of the channel matrices H2H2, H{Hi, and 
HlH! with Hi G CM nsXnd , respectively. 

Proof: See Appendix [Vj □ 

IV. Performance Analysis 

In this section, we use the new statistical results of A max derived in the previous section to investigate 
three important performance measures for beamforming schemes in MIMO dual-hop AF systems, i.e., 
outage probability, average SER and ergodic capacity. 

A. Outage Probability 

The outage probability is an important quality of service measure, defined as the probability that the 
instantaneous SNR falls below a pre-defined threshold 7th, and can be expressed as 



7th 



'_iy a {t-n s ) det h±X\ 



Pout(7th) = Pr(7 < 7th) = Pr A max < = = r — — -A- ^ (28) 

Note that (1281) presents the exact outage probability expression and allows for efficient evaluation of the 
outage probability of the optimal beamforming MIMO dual-hop AF systems for arbitrary numbers of 
antennas at the source, relay and destination terminals, at any desired SNR. 

Fig. [2] investigates the joint impact of n s , n r , and on the system outage probability. Results indicate 
that while it is always beneficial to have more antennas at the transmitter than the receiver, the advantage 
of doing so heavily relies on n r . For small n r , e.g., n r = 1, the benefit of putting more antennas at the 
transmitter is considerable. However, when n r becomes large, e.g., n r = 10, such advantage is significantly 
reduced, and n s and n d play a rather symmetric role. 



We now examine in more detail the outage probability for the special cases where the number of 
antennas at one of the nodes is one. We first consider the case n s = 1, i.e., when the number of antennas 
at source node is one. From (1281 ) and Corollary [1] the outage probability is 

Pout(Tth) = V ni r , ^Tvrv 7 jTT\ " (29) 

To gain more insights, it is of great interest to look into the low outage regime (or the high SNR 
regime). The outage performance in this regime can be generally characterized by two key parameters, 
i.e., diversity order and array gain. The diversity order is defined as the slope of the outage probability 
curve versus SNR curve (plotted on log-log scale), as SNR grows large. It is worth pointing out in the 
relaying system that high SNR consideration implies that both the source and relay power gains increase 
simultaneously^ Therefore, we are interested in the scenario where both p and a increase to infinity with 
a fixed ratio, i.e., a = kp and p — > oo. 

At high SNR, from Theorem |3] the outage probability can be approximated as 

p r det(tt) /7th V nm 

^i*)~ r{q + 1)uliT ^_ i + mq _ i+1) { ap ) . W 

which shows that the system achieves diversity order of q. 

Now, we consider the case when q = 1, i.e., the number of antennas at relay or destination node is 
one. From (1281 ) and Corollary |2] the outage probability of the system can be expressed as 

^)--wg^0 o, fe)^^( 2 ^)- (31) 

Note, when n s — 1 and p — 1, Eq. (l3TT) reduces to the results presented in [19, Eq. (9)]. At high SNR, 
from Theorem @] the outage probability can be approximated as 



f <-E^o(T)(£yr( P -n s +i) 

mk m T(p)T(n s ) \ p 



thJ 



m 

7th 



WrTthJ / 7th 



if p > n s 
if p = n s 
if p <n s 



mk m T(p)T(n s ) \ p 

n™T{n s -p) fj^ 
mk m T(p)T{n s ) ^ p 

which shows that diversity order of m is achieved by the system. To see how the outage probability 
behaves when n s increases, we consider the case p < n s . When n s is increased by one, the outage 
probability is reduced by a factor of which indicates that increasing the number of antennas at the 

2 If either of the source or relay power gain is fixed, then the outage probability of the system reaches an error floor when the other power 
gain increases. Therefore, the diversity order achieved in such scenario is zero. 
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source is always beneficial, but the benefit diminishes as n s becomes large. 

Fig. |3] shows the outage probability of the system with different rid when n s = 3, n r = 1. We observe 
that the Monte Carlo simulation results are in prefect agreement with our exact analytical results. Also, we 
can see that the high SNR approximation results are quite accurate even at moderate SNRs, i.e., p = 20dB. 
Moreover, the outage curves suggest that diversity order of 2 is achieved for n d = 2, and 3 is achieved 
for n d = 3,4, as expected. 

B. SER 

In addition to the outage probability, average SER is also another important metric for characterizing the 
performance of a communication system. For most modulation formats, the average SER can be evaluated 
as [26] 



SER = E 7 



a x Q (a/2027) 



(32) 



where Q(-) is the Gaussian Q-function, and a\, a<i are modulation-specific constants. 
Using (1321) and the pdf of A max in Theorem |2l the average SER can be evaluated as 

™ . , (-l)"-('-"'> EL g - s+ i Jo" aiQ W^px) det ($,(*)) dx 

SERbf p = rvr~F( • i 1YTV ~VT\ • (33) 

With arbitrary numbers of antennas at the source, relay and destination terminals, it does not appear that 
the integral in (|33l) can be expressed in closed-form. However, numerical integration can be performed to 
calculate SER much more efficiently than is possible via Monte Carlo simulations. 

For the special cases where the number of antennas at one of the nodes is large, accurate approximations 
can be obtained using Theorem [7] and the results in [27]. Define S(ci, c 2 , pj^ as the average SER of optimal 
beamforming in uncorrected Rayleigh fading channels with c\ transmit antennas and c 2 receive antennas 
at average SNR p. Then, we have the following result. 

Proposition 1: The average SER of optimal beamforming systems in MIMO dual-hop channels can be 
approximated as 

{S (nd, n s , nr . : p ) if n r — > oo, 

Proposition [T] offers an efficient means to evaluate the average SER of the system, and as illustrated 
in Fig. HI the asymptotic approximations are very accurate even at reasonably small number of antennas, 
i.e., n r or rid = 10, and it performs particularly good in the low SNR regime, where it is observed that 

3 The closed-form expression of <S(ci, C2, p) is given in [21]. To avoid cumbersome notations, we do not give the explicit expression here. 
It is also worth noting that the average SNR here is slightly different to that in [27], and it is given by p — ap. 
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the approximation curves and Monte Carlo simulation curves almost overlap. Moreover, the numerical 
results also suggest that for a given number of transmit antennas, it is beneficial to deploy more antennas 
at the receiver than at the relay. 

It is well known that optimal beamforming scheme achieves full diversity of MIMO systems, i.e., 
diversity of Cic 2 is achieved for MIMO systems with c\ transmit antennas and c 2 receive antennas. Hence, 
another important implication of Theorem |7] and Proposition \T\ is that full diversity order of nsn ^ nd [$ 
achieved by optimal beamforming in MIMO dual-hop AF systems when the number of antennas at one 
of the nodes becomes large. Now, we consider the special cases, where one of the nodes has only one 
antenna, and derive exact closed-form expressions for the average SER. 

Proposition 2: When n s = 1, the average SER of beamforming schemes in MIMO dual-hop AF systems 
can be expressed as 

ser bf( p) - -^(-D- 



^nuni-i+mp-i+i) 




Proof: The average SER of the system can be alternatively computed via 




aVP [°° e~P x 
Jo V% 



SERbf(p) = 77^= / —f=-^(x)dx. (36) 



Substituting the cdf expression presented in Corollary [TJ the desired result can be obtained with the help 
of the integration relationship Eq. (1251 ). and the fact that / °° ^j=-dx = □ 
To derive the diversity order and array gain of optimal beamforming in MIMO dual-hop systems, we 
now analyze the SER performance in the high SNR regime. To this end, with the aid of Theorem HI we 
can invoke a general parameterized SER result from [28] and perform some basic algebraic manipulations 
to obtain the following high SNR SER expression: 

™ am?r(n a + l/2) det(tt) q 

SERbf(p) ~ 2 ^ {a2k)qT{q + 1} ^ r(p _ % + 1)T{q _ l + l) P . C37) 

which indicates that the system achieves diversity order of q, as expected. 

Proposition 3: When q — 1, the average SER of beamforming schemes in MIMO dual-hop AF systems 
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can be expressed as 

" (ap) k (a 2 + l/p) k + 1 / 



a 



2 



x ^ + I, 1 + ,_ p _ ; ,_i_). (38 ) 
Proof: The proof is straightforward and omitted. □ 
Similarly, the average SER at high SNR can be approximated as 



SERbf(p) 



( am?' r(n,+l/2) (7)a < r(p-n a +Q _„ s 

20F(fca2)™=r(p)r(n s +i) Z 5 u p > n s 

ain?r(p+l/2) ln(|£ 



rp- p if p = n Sl (39) 



a 1 n?r( P +l/2)r(n B - P ) p tf p < ^ 



I 20F(fca 2 )Pr(p+l)r(n s )' 

ainPr(p+l/2)r(n s -p) 
I 2^(ka 2 )PV{p+l)V{n s )l 

which illustrates that diversity order of m is achieved. It is also of interest to compare the SER performance 
between optimal beamforming and OSTBC. In [14], the high SNR approximation of SER was derived 
for the OSTBC system employing BPSK modulation as 

f (^)" a r(n a +|)sr = f (7)a < r( P -n a +o(^)-' 

SER stbc(p)= ^Ti^n i ^ 

I 2v^fcpr( P +i)r(n s ) P n p < n s , 

where i? denotes the code rate of OSTBC schemes. Comparing (|39l and (1401) . it can be seen that both 
optimal beamforming and OSTBC achieve the same diversity order m. Since beamforming scheme exploits 
channel knowledge at the transmitter while OSTBC dose not, one would therefore expect beamforming 
scheme to outperform OSTBC. To see this, let us assume SER B f(pbf) = SERostbc(postbc)> and quantify 
the advantage of beamforming scheme in terms of SNRs. For the cases p ^ n s , it is not difficult to observe 
that 

A d B = (POSTBc)dB - (PBF)dB = (-R^s)dB, (41) 

which indicates a significant power saving of beamforming schemes over OSTBC schemes. 

Fig. [5] compares the average SER performance of optimal beamforming and OSTBC. For the OSTBC 
simulations, a full-rate Alamouti code [2] was used, so that R = 1, and the analytical result was plotted as 
per [14, (26)] J3 Clearly, we see that optimal beamforming and OSTBC achieve the same diversity order of 
2. Moreover, result indicate that the optimal beamforming scheme outperforms OSTBC by 3dB in terms 
of SNR, as is expected. 

4 There is a typo in the expression in [14]: the second argument of the Hypergeometric U-function should be "p + I — k + 1" instead of 

"p + I - k - 1". 
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C. Ergodic Capacity 

In this section, we study the ergodic capacity achieved by the optimal beamforming scheme. The exact 
ergodic capacity of the system can be expressed as 

C{p) = \^ {log 2 (l + a P A max )} . (42) 
Using (1421) and the pdf of A max in Theorem |2l the average SER can be evaluated by 



n) 2nLir(<?-i + i)r(p-< + i) 



With arbitrary numbers of antennas at the source, relay and destination terminals, the integral in <143T > is 
not known to be evaluated in closed-form. However, numerical integration can be performed to calculate 
the ergodic capacity much more efficiently than is possible via Monte Carlo simulations. For the special 
cases where the number of antennas at one of the nodes becomes large, very accurate approximations can 
be obtained. 



Let X(ci, c 2 , pf\ denote the exact ergodic capacity achieved by optimal beamforming in uncorrected 
Rayleigh fading channels with c\ transmit antennas and c 2 receive antennas at average SNR p. Then, we 
have the following result. 

Proposition 4: The ergodic capacity achieved by optimal beamforming in MEMO dual-hop channels 
can be approximated as 

~ (X (n r , n d , n s p + a) — X (n r , n^, a)) as n s — > oo, 
C(p) w I \X (n d , n s , ^iP) as n r oo, (44) 

|X (n r , n s , -^^ip) as n rf ^ oo. 

Proof: The results for the last two cases can be obtained via straightforward application of Theorem 

[71 For the first case, the key observation is that the ergodic capacity of the system can be approximated 

as 

C(p) = {log 2 (1 + (n s p + a)A max i)} - {log 2 (1 + aA max i)} . (45) 

Then, the desired result follows immediately. □ 
Fig. [6] illustrates the ergodic capacity of optimal beamforming MIMO dual-hop AF systems. As we 
can see, the capacity approximations in Proposition |4] are accurate over the entire SNR range of interest 
with even moderate number of antennas. Results also show that employing more antennas at the source 

5 The closed-form expression of I(ci,C2,p) is given in [27]. 
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results in better capacity performance than at the receiver, which agrees with the observations made from 
the outage probability and SER analysis. 

The above proposition deals with the case where one of the nodes is equipped with a large number of 
antennas. Next, we consider the case where one of the nodes is equipped with only one antenna. When 
n s = 1 or q = 1, we can obtain a tight upper bound for the ergodic capacity via Jensen's inequality, i.e., 

C{p) <Cu P = \ log 2 (l + apE{\ max }). (46) 
Alternatively, the following simple and accurate capacity approximation can be used [29] 

C(p) « G, K = i.o g2 e (.„(! + a „E {W ) - «Vg{^-^W') ) m 

To obtain further insights, we look into the high SNR regimes, and derive simple and closed-form 
expressions for the two key parameters dictating the ergodic capacity, i.e., high SNR slope and power 
offset Coo. The high SNR slope Soo and the power offset can be obtained via 

Soo = lim ^L, = lim ( log 2 p - ^M) , (48) 

p^oo log 2 p P^oo y Ooo / 

respectively. 

From (|42l) . we can easily establish that = \. Hence, the power offset can be computed as (the 
proof is given in Appendix IVTI) 

£oo,bf = log 2 j - log 2 e U(n s ) + V(p) - ^ e^E l+1 (^) J . (49) 

Note from the above expression that while the impact of p on is not immediately clear, increasing n s 
reduces C^, and hence improves the ergodic capacity of the system. More precisely, when the number 
of transmit antennas is increased by k\, £oq,bf will be reduced by log 2 e(ip(n s + kx) — ip(n s )), which is 
significant for small n s , but becomes almost negligible when n s is large. 

Fig. |7] illustrates the ergodic capacity of the system with difference antenna configurations. We observe 
that the capacity upper bound C up performs reasonably well over a wide range of SNR, especially at 
low SNRs, where C up almost coincides with the exact capacity. In addition, the high SNR approximation 
also provides a good reference to the exact result at even moderate SNRs, i.e., 5dB. Moreover, as n s 
increases, the performance of both C up and high SNR improve. For instance, the upper bound C up , high 
SNR approximation and the exact result C almost overlap when n s = 10. Furthermore, we see that the 
capacity approximation C 3pp is incredibly accurate for the entire SNR range, regardless of the numbers 



15 



of transmit antennas. 

V. Conclusion 

This paper examined the analytical performance of the optimal beamforming scheme in MIMO dual- 
hop AF systems. After deriving a host of new statistical properties of the maximum eigenvalue, exact 
and approximate closed-form expressions were derived for three important performance metrics: outage 
probability, average SER and ergodic capacity, which can be efficiently evaluated. Simple and informative 
expressions were also obtained for some special cases at high SNR, which illustrate the diversity order 
and array gains achieved by optimal beamforming. Our results have indicated the coupling effect of the 
numbers of antennas on the system performance. Moreover, optimal beamforming was shown to achieve 
the maximum diversity order of MIMO dual-hop AF systems as OSTBC schemes while providing a 
significant power gain. 
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Appendix I 
Proof of Theorem Q] 

We first present the following lemma, which will be used in the main proof of the theorem. 
Lemma 1: The statistical property of the maximum eigenvalue of H{h]> ^aH^Hg + I„ d j H 2 Hx is 
the same as that of the maximum eigenvalue of H|LHi, where Hi ~ CAf gi n 3 (0, l q g) I n J and L = 

A 2 

diag{/5 1? . . . , f3 q } with = , 2 ' + in which Aj is the singular value of H 2 . 

Proof: We first apply singular value decomposition on H 2 . Then, using the unitary invariant property 
of Gaussian random variables and [30, Lemma 1], we get the desired result. □ 

Now, we find it useful to consider two separate cases, namely: n s > q and n s < q. 

When n s > q, with the help of Lemma [H the cdf of the maximum eigenvalue A max conditioned on L 
was given in [31] as 

_ . ._. deb(^/-,(x)) 
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where Vi is a q x q matrix, with the determinant of 



and ^i(x) is a g x q matrix with entries given by 



= /3f- i+i 7 [n s - i + 1, , (52) 

where 7(72, x) is the lower incomplete gamma function. To obtain the unconditional cdf of A max , we must 
further average (|50T ) over the joint pdf of {0 < < ■ ■ ■ < f3 q < ^} [12] 

1 ^ a ^0Te-^ 

i<3 

Therefore, the unconditional cdf can be obtained as 



f(Pi, ...,fi 9 ) = fvr-r { ttttw ^TTT 11^ - II 77 (53) 



[*i(x)]ij= / /3^ +1 7 n s -i + 1,-)F->IL- -——dp. (56) 



n 



^A max (a:) = E L { 7* Amax (z I L) } = ; , 1 — , — — , (54) 

1 Ii =1 r(n s - % + i)r (g - % + 1) r (p - % + 1) 

where 

2i=/ detOM*)) (A-^-)^j ^t— (55) 

h<H<-<fi q <\ l<i<j<q [i ~ aPi) 

Utilizing the technique proposed in [32], (|55l ) can be evaluated as X\ = det(<&i(x)), where <&i(x) is a 
5x5 matrix with entries defined as 

3 

-a, 

(3J^ (1 - a/3)p+^ 
Making a change of variable y = y^g, the above integral can be computed as 

[<&i(x)k, = y p+ ^(l + ayf^e-^ (n s - i + 1, ^1±M^ (57 ) 

To this end, utilizing the series representation of incomplete gamma function [25, Eq. (8.351.1)], we have 

[$i(x)ki = TK - i + l)(4i - e-*Bi(i)) ( (58) 

where 

Ax= / ^- i -^(l + at) i+J '- 2 e-*tZt= V . )a l T(p + q + 1 - i - j + 1), (59) 
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and 



s ™fe poo 

B ^ x ) = J2lJ t^-^-^l + at) i+j+k - 2 e- t -*dt 
k=o k - Jo 

n s -i u i+j+k-2 /...., \ 

>r^X K X - (l + J+ k-2\ l P+q +i-i-.i-k+l _ 

= 2 Z^7J" 2^ I / ) aa; 2 ^ P+ ,+i-i-i-fc+i(2 v ^), (60) 

where the final expression is obtained using the integration relationship [25, Eq. (3.471.9)]. 

Now, we proceed to consider the case n s < q. The cdf of the maximum eigenvalue of F conditioned 
on L is given by [31] 



'_ 1 )n.(9-n.) det(* 2 (z)) 
dettVOn^rK-i + l) 



^a_(^|l) = 3 , rtr \ nn , r ,; ^ :^ (6i) 



where SE^rc) is a g x g matrix with entries given by 



\ q-n s -i 

— for i < q — n s , 

[*2(x)]ij={ V , , (62) 

Pf l+ \(q-i + l,f:) for i>q-n s . 

Hence, the unconditional cdf can be obtained by 

^ JX) = El ^^" L )> = m i T(n s - i + l)Ul l T(q- i + l)T( P - l + iy (63) 



where 

Jo<^ <...<fi q < i u - w 

which can be computed as X 2 = det($ 2 (x)), where <fr 2 (x) is a g x g matrix with entries defined as 



(64) 



f . , , loM-l) g — PP+ [Z:^ W for i <q-n s , 

J? /^ +1 7 (? - < + 1, f ) for . > g - n.. 

Following similar lines as in the first case, the integrals can be explicitly solved and we have 



^y-ns-i Y,^' 1 { q ' l+ l 3 - l )a l Y(p + i + j) for i<q- 



[* 2 (*)k,- = < ' ^' =u v ' ' (66) 

T(g - i + 1)(A 2 - e~ ax B 2 {x)) for i > g - n s . 

where 

A 2 = £ f s + ? + f ^ )aT(2 g + p + l-»-j-n fl + l), (67) 
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and 

9~* k 
k=0 

n s +i+j+k-q-2 , , • , • , i o\ 

As such, the unified expression can be obtained by appropriately choosing the parameters. 

Appendix II 
Proof of Theorem [2] 



Noting that f\ max (x) = cLFx ^ x> , an d making using of Theorem Q] and a classical formula for the 
derivative of a determinant, the pdf of the maximum eigenvalue can be written as 

where $>i(x) is a q x g matrix with (z, j ) -th entry given by 

, \&(x)]ij for i ^ I, 
= < (70) 

ax 

The derivative can be computed as 



doc V d/0c 

To proceed, we compute the first derivative of with respect to x explicitly as 

dB(x) _^l T( ^ k fr(i,j) + k\ ^(x fc ( v ^) e (-^)^- fc + 1 K e(M)+ ,_ fc+1 (2 v ^)) 

dx - Z 2^ k \ 2^ \ i ) a dx • W 

k=0 i=o v 7 

With the help of the differential property of Bessel-K function 

= (73) 

we get 

rf(x fc ( v ^) e ^) + ^ fc+1 ^ (M)+ ,_ fc+1 (2 v / i)) 
fix 

fl(j,j)+i-fc+i , efjjy+i— k 

= kx K L x 2 K e(i j )+l _ k+1 (2y/x) -x K x 2 ir e(i!j)+/ _ fc (2Vx). (74) 

As a result, some simple algebraic manipulations yield the desired result. 
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Appendix III 
Proof of Theorem [3] 

When n s = 1, A max is statistically equivalent to b^Ah, where h e CJ\f qxl and A = diag{/3i, . . . , /3 q } 



is a 9 x g diagonal matrix with for i = l,...,q being the g non-zero eigenvalues of the matrix 

Conditioned on A, A max is the sum of q independent chi-square random variables. Hence, the near zero 
behavior of A max can be expressed as [33] 

/A^A^^^n^v- 1 . (75) 

For this reason, the unconditional asymptotic expansion can be obtained by averaging over A. Applying 
the joint pdf of presented in [12], we have 

i<j i=l 

where the multiple integrals can be solved using the technique in [32]. Hence, 



m uu r(p -1+ i)r( g -i+i) L""L" n<ft- - n ^jS> "A ™ 



x 9 - 1 det(X) 

where X is a g x g matrix with entries 

pC]ij= [" P p - q+l+j - 3 e-^(l-a(3)- p ~ q d(3. (78) 
Explicitly evaluating the above integral yields the desired result. 

Appendix IV 
Proof of Theorem [4] 

To obtain the asymptotic first order expansion, we find it more convenient to start from the following 
alternative equivalent pdf expression of A max presented in [12]. 

which can be written as a Maclaurin series, and we are interested in the first term with the minimum 
exponent and non-zero coefficient. To achieve this, we consider three separate cases: 
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1) p > n s . We first notice that the minimum exponent for x ns ~ l e~ ax is n s — 1. Hence, the next key 
step is to determine the minimum exponent of the summation in (|79|) . When p > n s , it is easy to see that 
p — n s + i > 0, for i — 0, . . . , n s . Then, based on the asymptotic series expansion of K v (x), we have 

p-n s +i-l | - t \ 

x 2 ir p _ ns+J (2v^) = - ^ — r(fc + i) — ( ^ 



In x - V>(fc + 1) - tp(p - n s + i + k + 1) fc 



2 ^ r(A; + 1)I> - n s + i + A; + 1) 



Clearly, the dominant term is given by k = in the first term, i.e., 



x^K n+l (2^) « ^ S (81) 



Therefore, 

/u * , " f iS(^-" ,+i) ' <82) 

2) p = n s . In this case, the summation in (1791) can be expressed as 

^ I n . s Ja i xtK +i (2 v ^) = if (2v^) + ( U - W^Ki{2Vx). (83) 

i=0 ^ 1 ' i=l ^ 1 ' 

The dominant term can be found using (1801) and the following asymptotical expansion of K (x) 

^o(^)~-c-ln|, (84) 

where c is the Euler constant, and we therefore have 

x n s -i ( ns fn\ \ 

h - <*> * fww) p - ln 1 + S ( /) a ' r(p - n - + i] ) ■ <85) 

3) p < n s . Since the modified Bessel function of second kind is a symmetric function with respect to 
v, i.e., K v (x) = K- V (x), we can express the summation in (1791) as 

n a —p—l 



n B , v iis—p—i. , v 

W^Wa^iW^H £ ( ^ ) a^x^ K n ^ p . t (2^) + 

8=0 ^ ' 1=0 ^ ' 

f"* )a--^ (2v^)+ £ hy^^K^^). (86) 
A close observation shows that the term with minimum exponent comes from the first part in the right 
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hand side of the above equation. Hence, as x — > 0, we have 

Us—p—l 



4=0 ^ ' 4=0 ^ ^ 

« "gV-^fMa' ^-/" 1 ' 1 . (87) 
Clearly, the dominant term with minimum exponent is associated with i = 0. Hence, we have 

£ ^yx-^K p _ ns+l {2^) - I^^lx— . (88) 



Therefore, we get 



f^M - W^A^ 1 - (89) 

r(p)r(n a ) 



Appendix V 
Proof of Theorem [7] 

When n s — > oo, by law of large numbers, we have 



lit 1 1 = I. (90) 



n s -5-oo n* 



Noting that A max is the maximum eigenvalue of H 2 ^aH 2 H 2 + I Tld j H 2 HxH{, A max is equivalent to the 
maximum eigenvalue of n s H.\ ^aH 2 H 2 + l nd J H 2 , which completes the first claim. 
When rid — > oo, by law of large numbers, we have 



lim — ^ = I. (91) 



n d ->oo 

Then, we see that H 2 ^aH 2 H 2 + l nd j H 2 ps a nJ+i ^> anc ^ me secon d claim follows immediately. 
When n r — > oo, by law of large numbers, we have 



lim ^ = I, (92) 

-1 



and a close observation reveals that H 2 ^aH 2 H 2 + l n ^j H 2 has only rid non-zero eigenvalues which is 
identical to aw " r +1 . Thus, applying [30, Lemma 1] yields the desired result. 



22 



Appendix VI 
Proof of (1491) 

From the definition, we can compute the high SNR power offset as 

11 

£oo = log 2 y - E{log 2 A max }. (93) 
When q — 1, the maximum eigenvalue can be alternatively expressed as 

Ama , = (94) 

ah 2 h 2 + 1 

where hi G CA/" rasXl and h 2 G CAf pxl . Hence, we have 

E{log 2 A max } = E{log 2 (h;ha)} + E{log 2 (h 2 h 2 )} - E{log 2 (ah 2 h 2 + 1)}. (95) 
As such, the desired result can be obtained by using the results presented in [34,35]. 
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Fig. 1: The cdf results of A max : Analytical versus Monte Carlo simulations. 
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Fig. 2: Impact of n r on the outage probability of optimal beamforming MIMO dual-hop AF systems 

when a = p. 




Fig. 3: Outage probability results of optimal beamforming MIMO dual-hop AF systems when n s = 3, 

n r = 1, and a = 0.5p. 
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Fig. 4: Average SER of optimal beamforming MIMO dual-hop systems when a = p: Exact versus 

asymptotic approximations. 
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Fig. 5: Average SER of optimal beamforming MIMO dual-hop systems when n s = 2, n r = 1, = 3, 
and a = p, with BPSK modulation: Beamforming versus OSTBC. 
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Fig. 6: Ergodic capacity of optimal beamforming MIMO dual-hop AF systems when a = p: Exact 

versus asymptotic approximations. 
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Fig. 7: Ergodic capacity of optimal beamforming MIMO dual-hop AF systems for different n s , when 

n r = 1, rid = 4, and a = p. 



